Effects of finite superconducting coherence lengths and of phase gradients in 
topological SN and SNS junctions and rings 



D. Chevallier,! D. Sticlet/ P. Simon,i and C. Bena^'^'Q 

^Laboratoire de Physique des Solides, CNRS UMR-8502, 
Universite Paris Sud, 91405 Orsay Cedex, France 
'^Institute de Physique Theorique, CEA/Saclay, Orme des Merisiers, 91190 Gif-sur-Yvette Cedex, France 

(Dated: October 26, 2012) 

We study the effect of a finite proximity superconducting (SC) coherence length in SN and SNS 
junctions consisting of a semiconducting topological insulating wire whose ends are connected to 
either one or two s-wave superconductors. We find that such systems behave exactly as SN and SNS 
junctions made from a single wire for which some regions are sitting on top of superconductors, the 
size of the topological SC region being determined by the SC coherence length. We also analyze the 
effect of a non-perfect transmission at the NS interface on the spatial extension of the Majorana 
fermions. Moreover, we study the effects of continuous phase gradients in both an open and closed 
(ring) SNS junction. We find that such phase gradients play an important role in the spatial 
localization of the Majorana fermions. 

PACS numbers: 73.20.-r, 73.63.Nm, 74.45.+C, 74.50.+r, 



I. INTRODUCTION 

There has been a lot of recent interest in the possibil- 
ity to realize and detect Majorana fermionic states (see 
Refs. Handi for recent reviews). Following earlier the- 
oretical proposals semiconducting wires with a strong 
spin orbit coupling, such as InAs and InSb, in proximity 
of an s-wave superconductor have become a promising 
platform to observe Majorana fermions. The properties 
of such systems have been extensively studied for homo- 
geneous superconductor^ — but also in superconductor- 
normal (SN) junction^"— and superconductor- normal- 
superconductor (SNS) junctionsJ^ — 

A few recent experiments have reported the obser- 
vation of zero-bias conductance peaks the frac- 
tional ac Josephson effect in the form of doubled Shapiro 
steps^ both consistent with the existence of Majorana 
fermions. However such experimental observations have 
been shown to also have simpler explanations: a non- 
quantized zero bias peak could also arise in principle 
without any Majorana fermions^! — while the fractional 
ac Josephson effect can also occur in conventional SNS 
junctions and be understood simply in terms of Landau- 
Zener processes associated with the Andreev-bound-state 
spectrum of the junction. Therefore a non-equivocal 
proof for the observation of Majorana fermions is still 
lacking. 

Along the lines of finding the most appropriate ex- 
perimental strategy to detect such states, in a previous 
paper— we have analyzed the possibility to detect the 
Majorana states in NS junctions and SNS junctions us- 
ing local energy-resolved spectroscopy such as STM. We 
focus here on analyzing the effects of some relevant ex- 
perimental factors such as a non-zero penetration length, 
imperfect contacts, as well as a non-zero uniform phase 
gradient for the detection of such states in linear as well 
as ring Josephson junctions. 



In a previous work^ we have shown that the forma- 
tion of Majorana modes in a finite-size normal link can 
be controlled by varying the phase difference between two 
superconductors (SCs) in an SNS junction. Such Majo- 
rana modes stem from the intrinsic Andreev bound states 
residing in such link and thus may extend over the entire 
normal link. Such numerical prediction has been subse- 
quently confirmed by Klinovaja and Loss^^ who obtained 
a full analytical expression of the Majorana wave func- 
tion in an SN junction. These analyses have focused on 
sharp NS interfaces, for which the SC order parameter 
varies discontinuously at the interface. In the present 
work we consider first a non-zero penetration length of 
the SC in the normal state, and therefore a continuous 
decay of the SC order parameter close to the interface. 
We use a tight-binding exact-diagonalization technique 
to describe a semiconducting wire with strong spin-orbit 
coupling, end-contacted to s-wave SCs, in the presence 
of a Zeeman field^i^. We find that for small penetration 
lengths with respect to the length of the wire, a normal 
region is preserved, and the Majorana modes arising in 
this region have a significant spatial extension. However, 
for larger penetration lengths, the entire wire becomes 
a topological SC, and localized Majorana fermions form 
only at its ends. 

We also study the effect of a uniform gradient of the 
SC phase in both a linear and ring Josephson junction 
(see Fig. [T]). We find that, consistent with previous ob- 
servations on infinite wires, for small phase gradients the 
system remains fully topological, while for larger phase 
gradients, the portion of the system subject to the gra- 
dient enters into a gapless phase which does not support 
extra Majorana states besides the end Majorana states. 
For a Josephson ring we find that Majorana states can be 
induced in the normal region if the phase gradient is not 
too small, and if the total phase twist is close to an odd 
multiple of tt. The periodicity of the energy modes with 
the total phase twist has a small shift from 27r, shift that 
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goes to zero for large systems. Similar to the linear case, 
the system enters into a gapless phase when the phase 
gradient becomes larger than a critical value. 

The paper is organized as follows. In Sec. II we 
present the considered setups and the corresponding 
tight-binding models. In Sec. Ill we study the effect 
of the SC coherence length on the location and spatial 
extension of the Majorana fermions in the wire. Sec. IV 
discusses the influence of the transmission at the NS in- 
terface. Finally, in Sec. V, we focus on the effects of 
a phase gradient for linear and ring SNS junctions. We 
conclude in Sec. VI. 



II. MODEL 

In this section we present two models: one correspond- 
ing to linear SN and SNS junctions which will be used 
to understand the effect of the finite superconducting co- 
herence length and the effect of the transmission at the 
NS interface (Part A), and the second one corresponding 
to a SNS ring junction which will allow us to study the 
effect of a phase gradient (Part B). 



A. Linear NS and SNS junctions 

Our starting point is a semiconducting nanowire with 
strong Rashba spin-orbit (SO) coupling.— In the pres- 
ence of a Zeeman magnetic field and in the proximity of 
a SC, such a wire has been shown to exhibit end Ma- 
jorana fermionic states j^i^i^ We focus on using such a 
wire to make SN or SNS long junctions by end-contacting 
it to one or two s-wave SCs (See Fig. [T] a) and b)). In 
the semiconducting wire the superconductivity is induced 
only for distances from the NS interface smaller than the 
SC coherence length (CL); here the CL acts as a super- 
conducting penetration length (SPL), see Figs. [2] and [3] 
for an intuitive representation of the CL, as well as for a 
schematic depiction of the spatial dependence of the SC 
order parameter in the studied junctions. 

To be able to perform a numerical analysis of this sys- 
tem, we model it as a ID tight-binding chain^i^ 

^ 1 
H = ^c](t - /i + Vzaz)cj - -^Cj{t + iaay)cj^i + h.c. 

N 
N 

+ E^2(^)e(^-^)/«"e*'^^c,tc,, + h.c. 

(1) 

where the sum is performed only over the N sites be- 
longing to the semiconducting wire (we do not include 
the two s-wave SC sites in the model, but we assume 
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FIG. 1. Scheme of the proposed setups : a) a semiconduct- 
ing nanowire connected to an s-wave SC b) a semiconducting 
nanowire connected at both ends to s-wave SCs c) a semicon- 
ducting nanowire connected to an s-wave SC ring. 
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FIG. 2. The coherence length for a superconductor (green) 
corresponding to the "size" of a Cooper pair. This CL acts 
as a SC penetration length in the normal part (gray). 



that the sole effect of the SCs is to induce a SC order 
parameter in the wire which decays exponentially with 
the distance from the NS interface). Also, ay and cr^ are 
Pauli matrices which act in the spin space, /i is the chem- 
ical potential, Vz the Zeeman field, a the SO coupling, 
Ai/2 are the SC gaps of the two SCs, and ^sc is the coher- 
ence length of the superconductors (we assume that the 
two SCs have the same CL). According to Ref. [sl, the 
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FIG. 3. A schematic depiction of the spatial dependence of 
the SC order parameter in NS and SNS junctions. The red 
dotted hnes indicate the 'soft' effective NS boundaries at a 
distance ^sc from the interface. 



factors (1 — ^) and ^ before the superconducting terms 
are introduced in order to smooth the transition of the 
effective superconducting gap between the two supercon- 
ductors. We consider a perfect transmission of the junc- 
tion(s) throughout all the calculations except for Sec.HVl 
For the NS junction depicted in Fig. [T]a) we have taken 
A2 = and for the SNS junction depicted in Fig. [T]b) 
we have taken Ai = A2 = A. We work in units of t = 1, 
and the SPL is given in units of the lattice constant 
a which we also take to be equal to 1. 



B. SNS ring junction 

We consider using a topogical semiconducting wire to 
make an SNS ring junction (See Fig. [He)). The super- 
conductivity is induced via the proximity effect in one 
part of the ring. This ring SNS junction can also be 
modeled using a tight-binding chain 

^ I 
T-L = ^^Cj{t - fi^ Vz(Jz)cj - - [c]{t + ia(jy)cjj^i + h.c] 



N2 



j=Ni 



(2) 



where ^ corresponds to summing over all the sites 
around the ring. The topological SC section lies between 
sites A^i and N2. All the other sites correspond to the 
non-SC semiconducting nanowire. The phase gradient 
is modeled as a V(/) phase variation between two neigh- 
boring sites j and j + 1 in the SC region. All the other 
parameters are identical to the ones presented in the pre- 
vious section. 



C. LDOS and MP reminder 

The quantities that we use to describe the above setups 
are the local density of states (LDOS) and the local x- 
axis Major ana polarization (MP) which are given by the 
following expressions^ 

AN 

ftH = E E S{cv-Ei)\aiJ (3) 
^;(^) = E E -^(^ - Ei}2B^(a^,/3is), (4) 

i=l 6=t,i 

where (a*^, a^^, Pj^) are the components of the wave 
function on site j for the i^^ eigenstate of the system in 
the {cj^^Cj^^Cj^^Cji) basis. An exact digonalization of 
the Hamiltonians ([T]) and (|2]) allows us to evaluate these 
two quantities. We focus on the limit when the system 
is in the topological phase by choosing Vz = 0.4, A = 
0.3, /i = 0,and a = 0.2. A finite width 0.0002^^^/^ for 
the delta functions is introduced in the numerical evalu- 
ations. 



III. EFFECT OF THE SC COHERENCE 
LENGTH ON THE LOCALIZATION OF 
MAJORANA FERMIONS 

In this section, we study NS and SNS junctions made 
of a semiconducting nanowire connected to either one 
or two superconductors. The superconductivity in the 
nanowire is induced solely via the proximity effect (see 
Fig. [2] and [3]), with the SC order parameter decreasing 
exponentially with the distance from the SC interface 
(oc e~^). 

NS junction We first study an NS junction (see Fig. [1] 
(a)). In Fig. H and El we plot the LDOS and the MP 
of the nanowire as a function of energy and position for 
various values of the penetration length. 

In Fig. m we observe the formation of a zero energy 
mode close to the left end of the nanowire, in the vicinity 
of the NS interface, for values of £^sc which are not too 
small (for very small values of (_sc^ for example (_sc = 
10, precursor Majorana states form at the end of the 
wire, but these states are not fully Majorana polarized, 
and their energy is small but non-zero, and varies with 
the different parameters in the system). Also, Majorana 
states form at the right end of the wire, but their spatial 
extension depends on the value of the SPL. Thus, for 
small values of the SPL, the right-hand-side zero mode is 
extended over a large region of the wire (see the first two 
panels in Fig. |4]). With increasing the value of the SPL, 
this zero-energy state becomes more and more localized 
towards the end of the wire. 

To confirm that these states correspond to actual Ma- 
jorana fermionic states we plot the MP along the x-axis. 
Indeed, Fig. [5] shows that the left-end zero-energy state is 



4 




20 40 60 80 100 
Site 



20 40 60 80 100 "o 

Site 




40 60 80 100 
Site 



40 60 80 100 "o 

Site 



FIG. 4. LDOS for an NS junction as a function of energy and 
position for various values of ^sc- From left to right and up 
to down: ^sc=50, 100, 400, 600. 



a Majorana fermion with a total polarization of —1 (same 
as m Ref.i and |13| we obtain the total MP by integrat- 
ing over the entire peak) while the right-end extended 
Majorana carries a +1 polarization, consistent with the 
conservation of the total MP in the nanowire. 




tion in Ref.[T3land[T6l, and by assuming that the topolog- 
ical SC region induced in the wire is roughly given by the 
size of the SC coherence length. Moreover one can con- 
sider that an effective 'soft' NS boundary is created at a 
distance from the NS interface equal to the SC coherence 
length (see Fig. [3]). Thus, when the penetration length is 
much smaller than the nanowire length, the NS boundary 
is close to the left end of the wire, and the right Majorana 
mode extends over the large entire normal region. When 
the coherence length increases, the effective SC region 
becomes larger, the soft NS boundary is moving towards 
the right end of the wire, and the right-hand Majorana 
mode becomes more and more localized. When ^sc is 
larger than the length of the wire /, the entire wire is in 
a topological SC state, the Majorana states are localized 
at the two ends of the nanowire, and the extension of the 
Majorana wavefunction is constant and independent of 
the length of the nanowire {S^sc > I > £,m where ^rn is the 
"size" of the Majorana wavefunction). 

SNS junction We now turn to the analysis of SNS junc- 
tions (see Fig. [T]b)). In order to be as clear as possible, 
we present a plot of the LDOS for an SNS junctions with 
a zero phase difference, and a plot of the MP for an SNS 
junction with a phase difference of tt. Thus, for a zero 
phase difference, in Fig. [6l we show the LDOS as a func- 
tion of energy and position for various values of the SPL 
(modeled as before by an SC order parameter exponen- 
tially decreasing (oc e~^) with the distance from the 
two NS interfaces). 
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FIG. 5. MP for an NS junction as a function of energy and 
position for various values of ^sc- From left to right and up 
to down: ^sc= 50, 100, 400, 600. 

The spatial extension of the Majorana states can be 
understood in relation with the analysis of the NS junc- 



FIG. 6. LDOS for an SNS junction as a function of energy 
and position for = and various values of ^sc- From left to 
right, and up to down: ^sc=20, 40, 70, 100. 

We can see that, similar to the NS junction, effective 
NS 'soft' boundaries form inside the normal region for 
small values of the SPL {(^sc < 1/2). In this regime. 



5 



the situation is equivalent to the 'hard' SNS junction 
described in Ref. [l^. Indeed we note the formation of 
two zero-energy states at the ends of the wire, plus the 
formation of finite-energy quasi- Andreev bound states in 
the effectively-normal region. The two Majorana states 
at the ends of the wire have opposite polarization (±1), 
with no MP in the bulk of the wire (not shown). When 
we increase the SPL, these Majorana states become more 
and more localized towards the two ends of the nanowire 
and their spatial extension is constant and independent 
of the nanowire length (^^c > ^ ^m)- The normal 
region is becoming smaller and smaller, and the critical 
value of the SPL for which the entire wire becomes SC is 
given by half of the nanowire length (see Figs. [21 and [3]). 
Above this value no normal region forms inside the wire. 

The situation is more complex when the phase dif- 
ference is non-zero. As described in Ref. |l3|, in a long 
SNS junction with a non-zero phase difference, quasi- 
Majorana modes are being form inside the normal region, 
and these modes reach a zero-energy and have a full Ma- 
jorana character when the phase difference is equal to tt. 
Here we expect this behavior to hold for small values of 
the coherence length (_sc < and indeed, as plotted in 
Fig. [71 and consistent with Ref. UH, for ^ = tt we observe 
the formation of the zero-energy modes in the center of 
the wire; we have checked that these modes have a total 
polarization of —2, opposite to the two 'parallel' Majo- 
rana modes at the two ends of the wire which have a 
polarization of +1 each. 



ence of = tt the effective order parameter penetrating 
from the two ends of the wire will always cancel in the 
center of the wire^ (we assume that the two SC order 
parameters and SPLs are identical for the two SCs, if 
it is not the case, the point where the effective SC or- 
der parameter will cancel will be situated at a different 
position in the wire). Thus one will always form a non- 
topological region inside the wire in the region where the 
effective value of the SC order parameter is small, and 
consequently two Majorana states inside this region with 
polarization opposite to that of those at the ends of the 
wire. 

The coherence length, which is given roughly by the ex- 
tension of a Cooper pair, can greatly vary, from a few mi- 
crons in an usual superconductor (Al) to a few nanome- 
ters in the case of cuprates. While nanowires a few hun- 
dred nanometers in length are easy to obtain (typically 
the minimum size is around oc 500nm), we can expect 
that the transition ^sc > t>e achieved in typical 

setups. Moreover, the setup we propose here in which the 
proximity effect is achieved via the end contact with a su- 
perconductor, rather by placing the nanowire on top of a 
superconductor, may be easier to realize, and thus have 
technological advantages in the realization of topological 
SNS junctions. 
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FIG. 7. MP for a SNS junction as a function of energy and 
position for a phase difference of = tt and various values of 
Isc- From left to right, and up to down: ^sc=20, 40, 70, 100. 

For larger SPLs the central normal region becomes 
smaller and smaller, and for SPLs larger than the crit- 
ical value {^sc > V^) iio purely- normal region exists in 
the middle of the wire. However, due to the phase differ- 



IV. EFFECT OF THE BACKSCATTERING AT 
THE NS INTERFACE ON THE EXTENSION OF 
THE MAJORANA MODES 



In this section we study a 'traditional' topological NS 
junction for which the topological superconductivity is 
realized by placing a part of a semiconducting wire on 
top of a superconductor. It has been shown^^ that for 
such a junction, a Majorana fermion is localized at one 
end of the topological SC, while the second is extended 
over the entire normal part. Here we study the effect 
of the NS interface transmission on the spatial extension 
of the Majorana fermions. In Fig. [8] we plot the MP as 
a function of position and the tunneling coupling at the 
interface. 

When the tunnel coupling is zero, we recover the sit- 
uation of a topological SC completely disconnected from 
the normal region; indeed in this limit the right end Ma- 
jorana becomes totally localized in the SC part, at the 
interface. When we increase the coupling between the 
SC and the normal region, the right-end Majorana state 
extends more and more in the normal part, so that for 
a perfect coupling the Majorana mode is uniformly ex- 
tended over the entire normal part, as it has been shown 
also in Ref. and UH. 
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FIG. 8. Majorana polarization as a function of position and 
tunnel coupling at the interface. The 80 first sites correspond 
to the topological part and the last 20 to the normal one. 



V. THE JOSEPHSON EFFECT IN THE 
PRESENCE OF A UNIFORM PHASE GRADIENT 

In a recent paper ^ it was shown that super currents in 
the bulk of the superconductor could be used in principle 
to manipulate the Majorana fermions. More precisely, it 
has been shown that a constant spatial gradient in the 
phase of the superconducting parameter can drive the 
system from a topological phase supporting Majorana 
fermions to a trivial phase without zero-energy bound 
modes. Here we use the tight-binding model to investi- 
gate the Josephson effect in linear and ring geometries in 
the presence of a phase gradient. 

For the linear geometry, we consider a system similar 
to the long SNS junction explored in Sec. Ill A[ However, 
in this section we consider that the superconducting gap, 
the spin-orbit coupling as well as the Zeeman field are 
uniform over the entire wire, the only parameter depend- 
ing on position being the superconducting phase. A typ- 
ical phase spatial dependence is presented in Fig. [9l such 
that there are two outer regions of constant phase, with 
a central region uniformly experiencing a phase twist. 
We are interested to find if the uniform phase variation 
may give rise to a physics similar to the one observed in 
a traditional SNS junctionJ^ We also consider a semi- 
conducting ring with a "normal" region, and the rest 
experiencing a superconducting proximity effect and an 
uniform phase gradient. 



A. Constant phase gradient in a wire 

Before launching into the numerical analysis for the 
tight-binding form of the Hamiltonian in Eq. (|2]), we 
present an analytical study of the effect of the constant 
phase gradient on the topological index, that allows one 
to predict the existence of a topological phase transi- 
tion. It is not entirely surprising that the gradient of the 



superconducting parameter can make the system switch 
between a topological trivial and nontrivial phase. For a 
uniform gradient, this can be readily understood in the 
limit of an infinite wire. As shown in Ref. |39|, the phase 
of the pairing term can be gauged away, with the effect of 
adding a gradient-dependent correction to the canonical 
momentum and a renormalization of the hopping param- 
eter, as well as of the spin-orbit coupling. For the infinite 
wire, the condition to have a topological phase is given 
by^ 



t 



li-2t- -{V<i)f 



a2(V^)2 



V? + |A|2 \ < 0. 



(5) 



When the bandwidth t is larger than the other param- 
eters of the system, the second term of the product is 
always positive. Thus, for a zero chemical potential /i, 
the critical phase gradient is the exact lattice analogue 
of the continuum expression determined in Ref. [39|: 



(V0)c 



1/2 -> 1/2 



(6) 



A phase gradient has a Cooper pair breaking effect and 
can close the superconducting gap at the Fermi momen- 
tum. This leads to a second critical value for the phase 
gradient (V(/))^z above which the bulk gap closes, and the 
system enters a gapless regime. Its exact value is deter- 
mined by numerically studying the closing of the gap for 
an infinite system that experiences a uniform phase gra- 
dient V^. The momentum space Bogoliubov-de Gennes 
Hamiltonian is given by 

k 

H = {t — ii)tz — [tcos(/c) + a^m.{k)(jy\ cos{\/ (j)/ 2)) Tz 
— [tsm{k) — acos{k)(jy] sin(V(/)/2) + VzCTz — \^\^x 

(7) 

For our set of system parameters {Vz = 0.4, A = 0.3, a = 
0.2, /i = 0), the critical gradient for which the system 
enters into a gapless phase is {V(j))gi ^ 0.27, while the 
topological condition yields a critical phase gradient of 
(V0)c — 1.57. Thus, for the parameters we consider, the 
system enters first into a gapless metallic region before 
reaching a gapped trivial phase. 

We want to use this type of arguments to study a non- 
uniform wire with a variable-length central region that 
experiences the twist of the superconducting phase (we 
denote this region as gradient region - GR). When the 
phase gradient is zero, the system is in a topological phase 
with Majorana modes at the ends. The question that we 
are trying to address is if the central GR can be driven 
into a trivial gapped phase such that Majorana fermions 
can form also at the interface with the topological phases. 
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FIG. 9. Variation of the superconducting phase in a topolog- 
ical SC wire: the phase is set to zero in the left region, and 
a constant gradient is considered to be induced in the central 
region by bulk supercurrents such that the phase reaches a 
value of TT in the left region of the wire. 



We perform a numerical analysis for a system of 100 
sites. We start with a small value for the phase gradi- 
ent. In Fig. [To] the phase gradient over the 20 central 
sites is 7r/20, smaller than {V(j))gi value, and the GR re- 
mains topological. Hence there are only two Majorana 
fermionic states which form at the ends of the wire. Be- 
cause of the relative phase (/) = tt between the two outer 
regions, the Majorana fermions have parallel Majorana 
polarizations. We should note that the formation of two 
identical Majorana modes at the two ends of the wire 
does not raise a problem from the perspective of the 
conservation of the Majorana polarization, since, quite 
puzzling, the zero-energy MP is compensated by a con- 
tinuous distribution of MP in the higher-energy bands (if 
we integrate the MP over energy at any given site, the 
integral yields zero). 



Majorana Polarization Pm^ 




FIG. 10. The phase of the superconducting parameter is 
twisted by tt over the 20 central sites. The central re- 
gions remains topologically nontrivial and only two Majorana 
fermions with the same polarization form at the ends of the 



Let us start with a wire in the topological phase. The 
goal is to create new Majoranas by driving the central 
GR into a trivial gapped phase using only an increasing 
uniform phase gradient. We find numerically that this is 
not possible, and that the system is always driven first 
into a gapless phase iVcfygi < {V(j))c- This result is con- 
sistent with the phase diagrams in Ref. [39|. Therefore no 
Majorana fermions form at the interface between the GR 
and the exterior topological regions. 

A special situation arises however when the gradient 
region is constructed as a series of Josephson junctions 
with a phase increase of (2n + l)7r between two neigh- 
boring sites, with n being an integer. The GR must 
consist of an odd number of sites, to ensure a relative 
phase difference {2n' + l)7r between the left and right re- 
gions. Such a situation can for instance be encountered in 
arrays of topological superconducting islands which has 
been shown in Ref. [40] to provide a macroscopic realiza- 
tion of the Kitaev model^i^. 

In this situation, extended Majorana fermions can 
form in the GR. In Fig. [TT] we represent such a situa- 
tion with two extended Majorana fermions forming for 
a phase difference of ^ = IItt over an 11 site GR. Inte- 
grating the Majorana polarization over the central region 
yields a total value of two, showing that only two Majo- 
rana states form in this region. 



Majorana Polarization 




FIG. 11. Phase gain of IItt over 11 sites. Two Majorana 
fermions form in the GR having an opposite polarization with 
respect to the end modes. 



Note that the system is invariant under a change of 
27r in the phase gradient. Hence for a gradient region of 
N sites, there is a 2t:N periodicity in the total relative 
phase (j) between the left and right ends of the wire. In the 
special case when the phase gradient is distributed over 
an odd number of sites with Majorana fermions forming 
at = (2n + l)7r, the periodicity in the total phase is 
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FIG. 12. Evolution of the first two positive eigenvalues with 
the phase difference (j) for a 11-site GR. The constant zero- 
energy mode corresponds to an end Majorana fermionic state. 
The energy of the other mode evolves with and reaches 
zero only when there is a (2n + 1)ti phase difference between 
two neighboring sites (we can see that this happens here for 
n = and n = 1, corresponding to = llvr and respectively 
= 337r). 
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FIG. 13. Evolution of the lowest-energy modes with the phase 
difference (j). Note the shift of the 27r periodicity to a 27T-\-2S(j). 

the ring (see Fig. [T4|) . Same as in the previous analysis, 
the MP of these extended modes is compensated by the 
MP distributed continuously in the high-energy bands. 



B. Ring with a constant phase gradient 

Here we investigate the presence of Majorana fermions 
in a ring geometry divided in two regions of variable 
length. One of the regions is in a normal state (does 
not experience a SC effect), and the second one is SC 
with a SC phase that may be twisted. We investigate 
under what conditions Majorana fermions do form in the 
normal region. 

In the absence of a phase gradient, there no Majorana 
fermions form inside the normal region. Even though the 
SC region is in a topological phase, the would-be Majo- 
rana modes at the two ends communicate through the 
normal region. Therefore they split and become normal 
fermionic modes lying in the SC gap. We have checked 
that for a non-zero but imperfect coupling between the 
normal and the SC region, quasi-Major ana modes of fi- 
nite Majorana polarization appear at finite energy in the 
normal region. The extension of these modes in the nor- 
mal region is smaller and smaller, their polarization goes 
to one and their energy goes to zero when the coupling 
between the normal and the SC region goes to zero. 

In what follows we will consider a perfect coupling be- 
tween the normal and the SC state. While no Majo- 
rana states are formed when (/) = 0, for peculiar values 
of the phase difference accumulated over the SC region, 
and for phase gradients that are not too large, Majorana 
fermions can form in the normal region. This can be seen 
in Fig. [13] where we plot the low-energy eigenvalues as a 
function of the total phase difference. We note that when 
this phase difference is slightly larger than tt one of the 
modes reaches zero-energy. We have checked that this 
is an actual Majorana mode by plotting the Majorana 
polarization as a function of position and energy across 



Majorana Polarization Pm^ 




20 40 60 80 100 



FIG. 14. Majorana polarization as a function of energy and 
position for a ring with a total phase difference of 7r + ^(/). The 
system is periodic with the first and last site being equivalent. 
Extended Majorana modes form and are distributed over the 
entire normal region. 

We have checked that the Scj) difference between tt and 
the phase corresponding to the formation of the Majo- 
rana mode, and consequently the shift of the eigenmode 
periodicity from 27r to 27r + 2(50, decreases with the size 
of the system and goes to zero for an infinite GR. 

This 2(7r + S(j)) periodicity in the phase difference (j) is 
preserved for gradients smaller than a critical value which 
for the chosen values of the parameters in our system 
is {V(j))gi = 0.27, corresponding to a total phase differ- 
ence of 21.67r. As described above, larger gradients 
are predicted to drive the GR to a gapless phase. How- 
ever our numerical simulations indicate that, while larger 
gradient values do seem indeed to take the system into 
the gapless phase, Majorana fermions may still survive 
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for some peculiar gradient values; we do not understand 
what is the origin of this phenomenon. 

0.08 I ' i ' i ' ^ ' ^ ' i ' i ' i ' 1 

0.06 - - 
0^ 0.04 - 




0.02 - 



0.0 I ' ' ' ' ' — ^ ' — ^ ' ' ' ' ' ' ' 1 

10 20 30 40 50 60 70 80 

GR width (sites) 

FIG. 15. The energy of the first two positive eigenvalues for 
a decreasing GR length. The total phase twist is kept to 
= TT, while the number of sites is varied. Close to the 
critical value of the gradient (V(/))^z = 0.27 ~ 7r/11.6, the 
system passes into a metallic regime. Due to finite-size effects 
our numerical analysis recovers this transition at a number of 
sites of 9 instead of 11 — 12. 

Another way to illustrate the evolution of the physics 
of the system with the value of the phase gradient is to 
fix the total phase gain = tt, and study the evolution of 
the low-energy modes with the number of sites in the GR 
(a larger number of sites is equivalent to a smaller phase 
gradient). As it has been illustrated in Fig. [131 due to 
finite size effects, one has an infinitesimally-low energy 
mode; we note that the energy of this mode becomes 
smaller and smaller with increasing the number of sites. 
When the number of sites is reduced we expect to reach 
the critical value of the phase gradient that signals the 
passing of the system into the gapless phase. For an 
infinite system this is (V0)^/ 0.27, corresponding to 
7r/11.6. We would thus expect to see a phase transition 
when GR reaches the size of 11 — 12 sites. Indeed a 
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crossing of the bands and lifting of the low energy modes 
happens for a size of the GR of about 9 sites (see Fig.[T5]). 

VI. CONCLUSION 

We have used a tight-binding approach to calculate the 
energy spectrum, the local density of states and the Ma- 
jorana polarization in topological SC wires in the pres- 
ence of spin-orbit coupling and a Zeeman field. We have 
considered the effects of a non-zero penetration length 
of the SC in the normal state, as well as those of a uni- 
form phase gradient in both linear and ring Josephson 
junctions. We have found that the ratio between the SC 
penetration length and the length of the wire can con- 
trol the spatial extension of the Majorana modes in the 
normal section. 

We have also observed numerically the transition be- 
tween the topological and gapless state for the section of 
the wire subject to a significant phase gradient, in both 
linear and ring junctions. For the linear junction, the 
only situation in which we have recovered extra Majo- 
rana states in the region subject to the phase gradient 
is when the phase gradient is chosen such that the phase 
difference between two adjacent sites is an odd multiple of 
TT. For the Josephson ring we have found that Majorana 
states can be induced in the normal region if the phase 
gradient is not too small, and if the total phase twist is 
close to an odd multiple of tt. An interesting observation 
for systems subject to phase gradients is the formation of 
zero-energy Majorana modes whose full Majorana polar- 
ization is compensated by a continuous Majorana polar- 
ization distribution in the high-energy bands, such that 
the energy integral of the MP on-site is equal to zero. 
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